REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  the  burden,  to  the  Department  of  Defense,  Executive  Services  and  Communications  Directorate  (0704-0188),  Respondents  should  be  av'are 
that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB 
control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ORGANIZATION. 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

20-FEB-201 1 


4.  TITLE  AND  SUBTITLE 

COMPUTATIONAL  ELECTROMAGNETICS 


3.  DATES  COVERED  (From  -  To) 

15MAR  2008  to  30NOV  2010 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

FA9550-08-1-0148 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

OSCAR  BRUNO 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

CALIFORNIA  INSTITUTE  OF  TECHNOLOGY 
1200  E  CALIFORNIA  BLVD,  MC  217-50 
PASADENA,  CA  91125 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

AF  OFFICE  OF  SCIENTIFIC  RESEARCH 
875  NORTH  RANDOLPH  STREET  ROOM  31 12 
ARLINGTON  VA  22203 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

UNLIMITED-A 


10.  SPONSOR/MONITOR  S  ACRONYM(S) 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

AFRL-OSR-VA-TR-2 012-0 650 


14.  ABSTRACT 

The  objective  of  this  basic  research  effort  is  to  provide  the  conceptual  foundation  for  improved  EM  computational  codes  that  can 
predict  the  propagation  of  EM  waves  through  complex  media  (turbulent  atmosphere  and  penetrable  platforms)  as  well  as  wire 
antennas.  The  PI  and  his  team  will  undertake  the  theoretical  and  numerical  descriptions  of  well-conditioned  numerical 
implementations  (both  surface  integral  equations  and  volumetric  discretizations)  of  Maxwell's  equations  in  various  settings.  These 
include  open  and  closed  surfaces  (both  smooth  and  rough),  curved  wires,  and  volumetric  propagation  with  variations  in  the 
refractive  index.  The  AF  needs  to  have  accurate  and  timely  predictions  of  EM  characteristics  of  various  objects  and  media  and  the 
research  pursued  here  will  speak  to  those  needs. 


16.  SECURITY  CLASSIFICATION  OF:  | 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

17.  LIMITATION  OF 
ABSTRACT 


18.  NUMBER  19a.  NAME  OF  RESPONSIBLE  PERSON 

OF  OSCAR  BRUNO 

PAGES  - 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

y  626-395-4548 


Standard  Form  298  (Rev.  8/93) 
Prescribed  bv  ANSI  Std.  Z39.18 


Computational  Electromagnetics 
Grant  #FA9550-08- 1-0148 
Final  Report 
February  20  201 1 
Oscar  P.  Bruno 


OBJECTIVES 


*  Volumetric  propagation  and  scattering 


*  New  Integral  Equation  Formulations  and  Solvers 


*  Offshoot  methodologies 


REPORT 

The  overall  work  in  the  proposed  areas  for  the  duration  of  this 
effort,  through  November  2010,  is  detailed  in  what  follows.  In 
particular,  our  efforts  gave  rise  to  1)  Volumetric  propagation  and 
scattering  solvers  for  propagation  through  thousand  of  kilometers  of 
realistic  three-dimensional  atmospheres  (completed  in  2009);  2) 
Development  of  high-order  solvers,  which  do  not  suffer  from 
restrictive  CFL  conditions,  for  solution  of  time-dependent 
wave-propagation  problems  in  inhomogeneous  media;  3)  Electromagnetic 
integral  solvers  for  open  surfaces  in  three  dimensions,  with 
applicability  to  design  of  reflectarray  antennas  and  negative -index 
materials;  4)  Fast,  high-order  solvers  for  wire  antenna  problems;  and 
5)  Solvers  for  general  PDE  problems  which,  not  based  on  use  of  either 
finite  elements  or  finite  differences  use  the  continuation  method 
instead,  and  have  been  shown  to  lead  to  unconditionally  stable 
numerics  for  a  wide  range  of  realistic  PDE  problems. 


*  VOLUMETRIC  PROPAGATION  AND  SCATTERING 

-  PROPAGATION  OF  HIGH-FREQUENCY  NONLINEAR  WAVES  THROUGH  TISSUE: 
Our  approach  to  this  problem  is  based  on  the  FC-AD  algorithms 
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mentioned  to  some  extent  in  the  second  subsection  of  the  section 
"offshoot  methodologies"  below.  Although  the  full  Navier-Stokes 
equations  for  tissue  can  be  adequately  reduced  to  a  second-order 
PDE  for  a  single  scalar  unknown,  recent  numerical  experiments 
(in  a  collaboration  between  Caltech’s  postdoctoral  associate 
N.  Albin  and  OB)  have  shown  that,  for  a  variety  of  reasons,  the 
first-order  system  of  Navier-Stokes  equations  is  better  suited 
for  use  in  the  present  context.  Indeed,  first  order  systems  do 
not  require  use  of  up- winding  strategies  (since,  as  we  have 
shown,  even  without  use  of  up- winding  or  other  special 
treatments,  the  spectral  FC-AD  approach  does  not  give  rise  to 
instabilities  when  used  in  connection  with  the  first-order 
system).  Further,  the  use  of  the  Navier-Stokes  equations  retains 
all  the  nonlinear  terms  which  for  the  high  amplitudes  in  the 
focal  region  of  a  HIFU  source  have  the  potential  to  be 
non-negligible.  Work  in  this  area  over  the  last  year  has  lead  to 
a  highly  efficient  parallel  solver  capable  of  evaluating  linear 
and  nonlinear  wave  fields,  as  well  as  a  solver  for  propagation 
in  media  with  variable  index  of  refraction. 

-  ATMOSPHERIC  PROPAGATION:  As  part  of  this  effort  we  developed  a 
solver  capable  to  evaluate  propagation  of  waves  across  thousand 
of  kilometers  of  realistic  three-dimensional  atmospheres:  the 
best  previous  solvers  were  restricted  to  two-dimensional  (range 
and  height)  refractive  index  variations.  The  numerical  method  we 
introduced  in  this  context  is  applicable  to 
non- spherically -symmetric  atmospheric  refractive  index 
distributions  that  amount  to  sufficiently  small  perturbations 
from  spherically  symmetric  smooth  distributions;  as  shown 
through  a  variety  of  examples  drawn  from  actual  atmospheric 
data,  the  departures  from  spherical  symmetry  that  actually  occur 
in  the  Earth’s  atmosphere  fall  within  the  domain  of 
applicability  of  the  proposed  methodology.  One  of  the  main 
elements  of  our  algorithm  is  the  well  known  Rytov  approximation: 
we  express  a  given  atmospheric  refractive  index  distribution 
n(r)  as  a  sum  nO(r)  +  nl(r)  of  an  "unperturbed"  spherically 
symmetric  and  smooth  refractive  index  nO(r)  and  a  small 
"perturbation"  nl(r)  which  contains  the  three-dimensional 
variations  of  n(r).  The  character  of  the  problem  under 
consideration,  however,  is  such  that  even  its  solution  on  the 
basis  of  Rytov’s  method  gives  rise  to  extremely  high 
computational  costs.  We  thus  resort  to  (and  extend)  a 
high-frequency  localization  methodology  introduced  recently 
[Bruno,  Geuzaine,  Monro  and  Reitich,  (2004),  Bruno  and  Geuzaine, 
(2006)],  which  reduces  computational  costs  in  the  high  frequency 
regime  through  localization  around  the  sets  of  points  of 
stationary  phase  (which,  as  shown  in  our  paper  with  J.  Chaubell 
and  C.  Ao),  actually  coincide  with  light  rays.  In  conjunction 
with  this  strategy,  we  evaluate  the  necessary 
spheric  ally -symmetric  Green’s  function  by  means  of  geometrical 
optics — which  is  permissible  in  view  of  our  assumption  of 
smoothness  of  the  underlying  unperturbed  atmosphere.  The 
solutions  of  the  geometrical  optics  problems  arising  from  both 
localization  and  evaluation  of  the  spherically-symmetric  Green’s 
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function  are  produced  through  high- order  evaluation  of  integrals 
and  differential  operators,  and,  in  particular,  do  not  require 
use  of  numerical  ODE  solvers.  Interestingly,  unlike  all  other 
numerical  methods  applicable  to  this  problem,  the  accuracy  and 
computational  costs  of  the  proposed  approach  do  not  change  as 
frequencies  are  increased.  Thus,  the  proposed  methodology  will 
become  even  more  attractive  as  it  is  applied  to  expected  higher 
frequency  (microwave)  applications  [Kursinski  et  ah,  (2002), 
Gorbunov,  (2007)].  The  validity  of  the  Rytov  approximation 
itself  as  a  solver  for  problems  of  wave  propagation  within  the 
atmosphere  has  been  the  subject  of  several  discussions  ([Brown, 
(1966),  Brown  (1967),  Fried,  (1967),  Keller,  (1969)]);  roughly 
speaking,  an  application  of  the  theory  of  [Brown,  (1966)] 
indicates  an  agreement  with  our  contention  of  validity  for  a 
tropospheric  region  of  the  order  of  1000  km  in  horizontal 
dimensions.  While  Brown’s  theory  ([Brown,  (1966)])  is  based  on 
statistical  assumptions  that  may  be  difficult  to  verify,  the 
results  presented  in  our  work  show  unequivocally  that  under  the 
types  of  refractive  index  variations  present  in  the  upper 
atmosphere,  the  Rytov  approximation  produces  very  accurate 
solutions  for  domains  of  thousands  of  kilometers  in  range. 


*  NEW  INTEGRAL  EQUATION  FORMULATIONS  AND  SOLVERS 

-  NEW  INTEGRAL  EQUATIONS  FOR  OPEN  SURFACES.  As  one  of  the  central 
parts  of  the  present  effort  we  are  developing  a  novel  approach 

for  the  solution  of  problems  of  diffraction  by  open  surfaces  in 
three  dimensional  space.  This  methodology  relies  on  use  of 
weighted  versions  of  the  classical  first-kind  integral  operators 
for  the  open- surface  electromagnetic  scattering  problems  (for 
which  the  unknowns  do  not  tend  to  infinity),  together  with  a 
generalization  to  the  open-surface  case  of  the  well  known 
closed- surface  Calderon  formulae.  We  had  previously  established 
for  the  acoustic  case  that,  when  used  in  conjunction  with  a 
spectrally  accurate  formulation  and  Krylov-subspace  linear 
algebra  solvers  such  as  GMRES,  the  Calderon-like  formulation 
produces  results  of  high  accuracy  in  very  small  numbers  of 
iterations  for  both  the  Dirichlet  and  Neumann  problems,  and  for 
low  and  high  frequencies  alike. 

-  SOLVERS  APPLICABLE  TO  HIGHLY  CHALLENGING  GRATING  DIFFRACTION 
PROBLEMS.  In  collaboration  with  M.  Haslam  we  introduced  a  new 

algorithm  for  evaluation  of  electromagnetic  and  acoustic 
scattering  by  one -dimensional,  perfectly  conducting  smooth 
periodic  surfaces  z  =  f(x).  Our  algorithm  is  applicable  to  both 
TE  and  TM  polarizations;  in  either  case  the  problem  is 
formulated  in  terms  of  a  second-kind  singular  integral 
equation.  Our  high- order  algorithm  for  these  problems  is  based 
on  two  main  elements:  (i)  Expansion  of  the  unknown  integral 
density  in  a  Floquet  series,  and,  (ii)  Super-  algebraically 
accurate  treatment  of  integrals  involving  singular  kernels.  For 
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an  incident  field  with  C 1  smoothness,  our  algorithms  converge 
faster  than  any  integer  power  of  M  and  N,  where  N  denotes  the 
number  of  surface -density  unknowns,  and  M  denotes  the  number  of 
weights  used  for  integration  of  the  kernel  functions.  Progress 
in  this  area  over  the  last  year  has  given  rise  to  solvers  for 
non- graph  diffractive  surfaces  containing  geometric 
singularities  such  as  corners  and  edges_.  The  latter  work 
resulted  in  the  WCRM  publication  listed  below. 

-  CURVED  WIRE  ANTENNAS  SOLVERS.  We  developed  an  extension  of  the 
methodology  previously  put  forward  for  straight  wire  antennas 

that  applies  to  arbitrary  curved  wire  antennas.  We  developed  a 
general  solver  for  open-ended  general  curved  wire  antennas,  and 
we  demonstrated  the  ability  of  our  solvers  to  evaluate,  with 
high  accuracy,  scattering  by  curved-wire  problems  hundreds  of 
wavelengths  in  size  in  computing  times  of  the  order  of 
minutes.  Further,  we  demonstrated  a  capability  to  solve  problems 
of  radiation  by  antennas  composed  of  large  number  of  antenna 
elements. 

-  PARALLELIZATION  OF  SOLVERS  FOR  ACOUSTIC  AND  ELECTROMAGNETIC 
SCATTERING.  Highly  efficient  algorithms  were  produced  which  can 

solve  problems  of  scattering  in  parallel  clusters.  For  example, 
our  solvers  demonstrated  the  largest  direct  solution  ever 
obtained  of  a  problems  of  scattering  -  containing  a  scatterer  of 
a  diameter  of  400  wavelengths,  and  with  an  accuracy  of  four 
digits,  in  a  13  hour  of  computation  in  256  processors  at  2.7  GHz. 


*  OFFSHOOT  METHODOLOGIES 


-  CONTINUATION-BASED  "ALTERNATING  DIRECTION"  SOLVERS  FOR  GENERAL 
PDEs:  FC-AD  ALGORITHMS.  A  new  methodology  was  introduced  for  the 

numerical  solution  of  Partial  Differential  Equations  in  general 
spatial  domains.  The  methodology  is  based  on  the  use  of  the 
well-known  Alternating  Direction  Implicit  (ADI)  approach  of 
Peaceman  and  Rachford  in  conjunction  with  one -dimensional  and 
high-order  accurate  Fourier  representations  of  non-periodic 
data,  obtained  by  way  of  the  "continuation  method"  we  introduced 
recently  in  the  context  of  the  problem  of  surface 
representation.  These  methodologies  were  further  extended  to 
enable  treatment  of  problems  involving  Neumann  boundary 
conditions,  variable  refractive  indexes,  and  a  parallel  code  was 
produced  on  the  basis  of  a  new  method  for  "overlapping  domain 
decomposition". 

-  CONTINUATION-BASED  "EXPLICIT"  SOLVERS  FOR  GENERAL  PDEs  and 
HYBRIDS:  A  variety  of  methods  for  solution  of  Partial 

Differential  Equations  on  the  basis  of  the  Fourier  Continuation 
method  were  developed,  and  are  the  subject  of  significant 
efforts  at  this  time.  Such  solvers  were  applied  to  a  variety  of 
problems  and  configurations,  including  the  Maxwell  system, 
propagation  of  elastic  waves  in  solids,  convection  diffussion 
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problems  in  tissue,  etc. 

-  COMPUTATIONAL  BOUNDARY  CONDITIONS.  Many  approaches  have  been 
proposed  for  this  fundamental  problem  in  the  field  of  wave 
scattering.  All  of  them  fall  into  one  of  three  main  categories. 

The  first  class  of  methods  is  based  on  mathematical 
approximations  or  physical  heuristics.  These  boundary  conditions 
are  often  local  in  space  and  time,  therefore  easy  to  implement 
and  run  in  short  computing  times.  However,  these  approaches  give 
rise  to  spurious  refleections  at  the  artificial  boundary,  no 
matter  how  refined  the  discretization  is,  which  travel  back  into 
the  computational  domain  and  corrupt  the  solution.  A  second 
group  consists  of  accurate  and  convergent  methods.  However, 
these  formulations  are  usually  nonlocal  in  time  and  space,  thus 
harder  to  implement  and  often  more  expensive  than  the 
computation  of  the  interior  scheme  itself.  Finally,  there  are 
methods  which  are  accurate  and  fast.  These  approaches  are  often 
local  in  time,  and  the  nonlocality  in  space  is  confined  to  a 
closed  surface  rather  than  the  whole  computational  domain.  The 
drawback  of  these  approaches  lies  in  the  fact  that  the  outer 
boundary  must  be  taken  to  be  either  a  sphere,  a  plane,  or  a 
cylinder.  For  many  applications  of  interest,  this  may  require 
use  of  a  computational  domain  much  larger  than  actually  needed, 
which  leads  to  an  expensive  overall  numerical  scheme.  In  our 
work  for  this  problem  we  introduced  a  new  methodology  in  order 
to  compute  the  fields  at  the  artificial  boundary.  The  boundary 
condition  is  both  nonlocal  in  space  and  time,  but  the  nonlocal 
behavior  is  confined  to  a  finite  number  of  points  in  time  and  to 
a  surface  in  space.  Like  the  second  class  of  methods  described 
above,  the  proposed  algorithm  is  accurate  and  numerically 
convergent,  yet  its  computational  cost  is  less  than  the 
underlying  portion  of  the  volumetric  calculation.  And,  unlike 
the  third  category,  this  new  approach  allows  us  to  choose  the 
artificial  boundary  to  be  arbitrarily  close  to  the 
scatterer.  This  method  is  based  on  a  concept  of  "sparse 
equivalent  sources"  representations  (introduced  by  the  PI  in 
collaboration  with  L.  Kunyansky)  which  allows  a  highly  accurate 
and  fast  evaluation  of  the  boundary  condition  when  used  in 
combination  with  fast  Fourier  transfonns.  A  variety  of 
numerical  results  we  produced  for  three  dimensional  problems 
demonstrate  the  accuracy,  high- order  ocnvergence  and 
effectiveness  of  the  methodology. 

RELEVANT  PUBLICATIONS 


"Mathematical  analysis  of  a  fast,  high  order  Nystrom  method  for  BIEs 
solution  of  three-dimensional  scattering  problems",  O.  Bruno, 

V.  Dominguez  and  F.  Sayas,  Preprint,  2011 

"A  highly  accurate  integral  method  for  open- surface  scattering 
requiring  small  numbers  of  GMRES  iterations",  O.  Bruno  and  Stephane 
Lintner,  In  preparation. 
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Plenary  speaker  at  the  Isaac  Newton  Institute:  conference  on  highly 
oscillatory  problems,  12—17  September  2010. 

Plenary  Speaker,  ICOSAHOM  2009,  Norwegian  University  of  Science  and 
Technology  in  Trondheim,  Norway,  July  22-26,  2009 

Plenary  Speaker  at  the  Oberwolfach  meeting  on  "Analysis  of  boundary 
element  methods";  Martin  Costabel  and  Ernst  P.  Stephan,  Organizers, 

April  13th  -  19th,  2008. 

Plenary  Speaker,  Workshop  on  wave  scattering  problems,  France  , 
December  12-14,  2007 

Organizer  of  the  workshop  "High-order  methods  for  computational  wave 
propagation  and  scattering"  (w.  R.  Kress)  at  the  American  Institute  of 
Mathematics,  Palo  Alto,  CA,  2007. 

Plenary  Speaker,  Leading  lecture,  at  the  Isaac  Newton  Institute 
Workshop  Effective  Computational  Methods  for  Highly  Oscillatory 
Problems:  The  Interplay  between  Mathematical  Theory  and  Applications  2 
July  to  6  July  2007;  Isaac  Newton  Institute,  Cambridge  CB3  0EH, 

Plenary  Speaker:  Workshop  on  Oscillatory  integrals  and  integral 
equations  in  high  frequency  scattering  and  wave  propagation",  June 
19,  2007.  A  one-day  workshop  at  the  Isaac  Newton  Institute  20  Clarkson 
Road,  Cambridge  CB3  0EH,  UK. 

Plenary  Speaker  at  the  Oberwolfach  meeting  on  Computational 
Electromagnetism  and  Acoustics,  February  5  to  9th,  2007. 

Plenary  Speaker  at  the  Seventh  Mississippi  State  -  UAB  Conference  on 
Differential  Equations  \&  Computational  Simulations  November  1-3,  2007 


file:///Q|/SHARE/Rcports%20Processing/DTIC%20Reports%20For%20Yvorme/Melvin%20List/FA9550-08-l-0148.txt[8/9/2012  11:30:00  AM] 


Doubletree  Hotel  Birmingham,  AL,  USA. 
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